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FOREWORD 


This report was prepared by Eagle Engineering, Inc., Hampton 
Division, under contract NASW 4430, sponsored by NASA Langley 
Research Center and the Federal Aviation Administration Technical 
Center under the FAA-NASA Interagency Agreement No. DTFA03-890- 
A-00019 of 13 June 1989. This report fulfills the requirement 
of the Program Plan for the National Aging Aircraft Research 
Program, DOT/FAA/CT-88/32 , August 1989, Paragraph 2 . 3 . 2 . 1 , Flight 
Loads. 

The Eagle Engineering, Inc. effort was performed by Norman L. 
Crabill and administered under the direction of Joseph W. Stickle 
(NASA Langley Research Center) and Thomas DeFiore (FAA Technical 

Center) . 
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THE NASA DIGITAL VGH PROGRAM- 
EXPLORATION OF METHODS AND FINAL RESULTS 
Volume I: Development of Methods 

Normal L. Crabill 
Eagle Engineering, Inc. 

Hampton Division 

SUMMARY 

The results of a NASA effort to utilize data from existing 
Digital Flight Data Recorders on airline transport aircraft in 
routine airline service indicates that many statistical data types 
useful to aircraft designers and operators can be compiled from the 
limited measurement types selected. Techniques for solving the 
significant problem of data editing were developed, along with 
methods for separating maneuver and gust accelerations using 
200 hours of data taken from L 1011 operations in 1973. Some 
results indicate that the acceleration derived exceedances at the 
4 samples per second rate generally available from the Digital 
Flight Data Recorders may be only 1/2 to 1/3 those obtainable at 
20 to 40 samples per second and, therefore, the present accelera- 
tion data must be used with caution. These techniques and methods 
are described and the results are given in Volume I. Similar 
analysis techniques were applied to about 5000 total hours for 
L 1011, B 727, B 747 and DC 10 aircraft operations of 1978 through 
1982. These results are given in Volumes II, HI/ IV, and V, 


respectively . 


INTRODUCTION 


The NACA-NASA has long had an involvement in determining 
actual operating conditions of commercial aircraft to aid designers 
in developing satisfactory design criteria. Starting in 1933, the 
NACA VG Program, using the smoked glass and stylus technique 
(ref. 1) , gave an analog representation of the operating VG diagram 
for direct comparison with the designer's load factor assumptions. 
A new "VGH" recorder was introduced in 1946 to give time histories 
of velocity, "G" load, and height, which had to be manually 
manipulated into suitable statistical forms to provide meaningful 
guidance to designers (refs. 1, 2). These programs involved many 
different types of aircraft, including general aviation and airline 
transport aircraft (refs. 3, 4). After 1971, however, the NASA VGH 
program was restricted to general aviation operations only 
(ref. 5). In 1977, with the cooperation of the manufacturers and 
airline operators, the NASA renewed the Data Recording Program to 
study airline operations and investigate the utilization of the 
Digital Flight Data Recorders (DFDRs) (ref. 6) required on all 
inrge turbine aircraft certificated since September 30, 1969. 
Limited early results of this effort are given in reference 7 . 

This report provides the final results of this exploratory 
Digital VGH ( DVGH ) program. Parameters utilized were a subset of 
those already available on the existing DFDRs without imposing any 
new requirements on data quality. No new recording system was 
used. The data quality problems that were encountered were 
handled with appropriate editing techniques. 
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The first volume of this report describes these data editing 
techniques, the analysis methods, and the many statistical data 
types developed in consultation with the airframe manufacturing 
industry using 200 hours of DFDR data taken from routine airline 
operations of a Lockheed L 1011-1 aircraft in 1973. Similar 
analysis techniques were subsequently applied to about 5000 total 
hours obtained from the L 1011, Boeing's B 727 and B 747, and the 
Douglas DC 10 aircraft in airline operations from 1978 through 
1982. These results are given in Volumes II, III, IV, and V of the 
present paper. 

Starting in 1982, NASA developed and flight-tested on 
contract, a brassboard version of a "Smart" Flight Recorder in 
which statistical data were computed in near real-time and stored 
on-board in the recorder (ref. 8). This approach eliminated the 
tedious manual labor required to edit-out the data errors intro- 
duced by the frequent off-nominal performance of the transcription 
process used in the DFDR system. This Smart Flight Recorder 
approach shows promise for obtaining large quantities of statisti- 
cal data for this and similar applications. These techniques and 
the results obtained in 200 hours of operation on a Beechcraft King 
Air are not discussed further in this report since they are 
described in reference 8. 
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SYMBOLS 


a 


n 


a n +1.0 


a 


nG 


a 


nM 


a 


y 


B 

c 

CAS 

eg 

c i 

L a 

dB 

DC 

deg 

DFDR 


DVGH 


r: 

L 


FLP 

ft 

g 

GW 


incremental component of normal acceleration near the 
aircraft eg; g units positive toward the top of the 
aircraft 

total normal acceleration near the aircraft c.g.; 
g units; positive toward the top of the aircraft 

incremental normal acceleration, identified as due to 
gusts, g units; positive toward the top of the aircraft 

incremental normal acceleration, identified as due to 
maneuvers, g units; positive toward the top of the 
aircraft 

lateral acceleration near the aircraft eg; g units; 
positive toward the right wing tip. 

Boeing 

aircraft mean wing chord, feet 
calibrated airspeed, knots 
center of gravity 

aircraft lift-curve slope, per radian as used in the 
equations 

decibels 

Douglas Commercial 
degree 

Digital Flight Data Recorder 
digital VGH 

frequency, cycles per second 
Federal Aviation Regulations 
flap 
feet 

acceleration of gravity 
gross weight, pounds 
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HP 

HZ 

kft 

K 9 

klbs 

kts 

L 

LaRC 

M 

NACA 

NASA 

R 

S 

SPL 

U de 

VG 

VGH 

US 

V e 

P 

Po 

^9 


> 


pressure altitude, feet 
frequency in cycles per second 
thousands of feet 

= Q.88 u . gust alleviation factor from reference 13 
5.3 +/Jg 

thousands of pounds 
knots 

Lockheed or Left; usage is obvious 
Langley Research Center 
Mach number 

National Advisory Committee for Aeronautics 
National Aeronautics and Space Administration 
right 

aircraft reference wing area, square feet 
spoiler 

derived equivalent gust velocity, feet per second 
velocity and "G" load measurement system 
velocity, "G" load, and height measurement system 

United States 

equivalent airspeed = True Air Speed x P — 

. * p o 

air density, slugs per cubic feet 

standard atmosphere sea level air density, slugs per 
cubic feet 

2 W , mass parameter from reference 13 

C La P cgS 

absolute value 

greater than 
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AIRCRAFT AMD INSTRUMENTATION 

Aircraft 

DFDR data were purchased from a U.S. airline operating a 
Lockheed L 1011-1 with three Rolls Royce RB-211-22 high-bypass- 
ratio turbofan engines of 42,000-pounds thrust each. A three-view 
drawing of the aircraft showing dimensions and control locations 
is given in figure 1. Weight and geometric characteristics are 
listed in Table I and the untrimmed flexible airplane lift-curve 
slope data used in the analysis are given in Table II. 

Instrumentation 

No new instrumentation was added to the aircraft. Instead, 

it was decided to access a small subset of the data already being 

obtained on the existing DFDR required by the FAR 121.343 and 

described in reference 6. Parameters finally selected by NASA, 

after discussions with the aircraft manufacturers, are given below 

with their range, sample rate, and accuracy taken from reference 6. 

Samples 
per 

Parameter Range and Units Second Accuracy 

a n +1.0 -3g to +6g 4 +.2 g’s stabilized, 

+10% transient 
(see page 6) 

a y -lg to +lg 4 +.05 g's stabilized, 

+10% transient 
(see page 6) 

4 +10kts 

1 +100 to 700 ft 

1 +3° 

1 unknown 

1 unknown 

i 

i 


CAS 


100 to 450 kts 

HP 


-1000 to 50,000 

FLP 


-5° to 60° 

SPL 

2R 

-5° to 60° 

SPL 

5L 

-5° to 60° 
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Accelerometer transient response is specified in reference 6 as: 


"8.4 Filtering (Output Frequency Response) 

The accelerometer should contain effective filter- 
ing means to screen out undesired high frequency 
vibration data. The output signal level should be 
3 dB below the signal levels set for in Section 
8.3.1 for vibration having a frequency of 4 Hz . At 
higher frequencies, the output signal levels should 
continue to decrease at the rate of 12 dB per 
octave . " 

The signal levels referred to in Section 8.3.1 of reference 6 are 
the null signal levels. 

Appendix A shows that in one test, normal accelerations 
measured at 4 samples per second with a DFDR system, like that 
described in reference 6, correlated well with normal accelerations 
measured with a more accurate NASA data system. In this test, the 
standard deviation of the DFDR normal acceleration measurement with 
respect to that of the NASA system was 0.055 g's. Thus, the normal 
acceleration measurement from the DFDR can be expected to be well 
within the ±0.2g quoted at 4 samples per second. Appendix A also 
shows, however, some unpublished results of a NASA B 57B flight 
test which indicate that for that test where the accelerometer cut- 
off was 10 Hz at 4 samples per second, the acceleration exceedances 
are 1/2 to 1/3 those obtained at 20 to 40 samples per second. 
Similar considerations may apply to the lateral acceleration 
exceedances. This possible limitation on the acceleration data 
presented herein, and in the companion Volumes II, III, IV, and V, 
must be considered in the application of the present results. 
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SCOPE OF DATA 


The data were obtained from flight operations of one regularly 
scheduled airline operating over the route structure from February 
1973 to May 1973, shown in figure 2. Some gaps in coverage did 
occur due to the characteristics of the DFDR "Crash Recorder" 
system which acts as a 2 5-hour loop tape erasing any data older 
than 2 5 hours. The aircraft flew 8 to 12 hours/day; thus, the 
airline company was asked to provide data from the recorder every 
2 or 3 days to provide continuity of data on one aircraft. This 
was not always achieved, however, and gaps in continuity did 
appear. Furthermore, some whole tapes were rejected in the edit 
process, resulting in additional gaps in the coverage of the 
service record. This was not considered serious in the development 
of this prototype data system. 

Eighty-three flights were utilized for a total of about 
200 hours and 91,000 nautical miles. The data were obtained by 
the airline operator in 1973 and were purchased by NASA in 1977. 
Subsequent development of the reduction techniques, including 
editing and definition of the final statistical formats shown 
herein, involved several interactions with airliner manufacturers 
and operators. 
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DATA REDUCTION 


Process 

The basic data reduction process is shown in figures 3, 4, 
and 5. The airline company supplied transcription tapes of the 
requested parameters, takeoff and landing gross and fuel weights, 
strip charts, and listings of the first 300 seconds of all of the 
parameters. The transcription tapes were converted into NASA 
engineering units tapes from which time history plots were made for 
editing purposes. 


Editing 

Originally, the time histories were "raw" or unedited and the 
editing was done manually. Later, a computer algorithm was 

developed to assist, and after some experience was gained in its 
use, the edit program was used to produce the engineering units 
time history plots. Appendix B describes the development of the 
KEDIT program and its application. It is a two parameter local 
signature analysis program designed to "pluck" or remove wild 
joints and replace them with reasonable values — it is not a 
-moother. Other edit functions manually performed include: 

1. Identification of lift-off and touchdown times 

2. Bias removal 

3. Data overlap deletions 

4. KEDIT performance 

5. Acceptance or rejection of each flight 

An example of the results of this process is given in figure 6. 
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Weight Calculation 


The weight at any time is found by a linear interpolation 
between the gross weight at lift-off and touchdown. Calculations 
summarized in Appendix C using manufacturer's fuel-flow equations 
indicated that the discrepancy was, at most, 2 percent occurring 
at the top of the initial climb on long flights. The difference 
was considered small and not worth the effort to reduce it. 

Counting Technique for Accelerations 
It was decided to utilize the level-crossing counting 
technique for acceleration exceedance analysis, due to its ease of 
application on the computer and its preference for design criteria 
evolution as opposed to fatigue life tracking, as pointed out in 
reference 9. The previous analog VGH programs (refs. 1, 2, 3, 4), 
using manual and eyeball methods, of necessity employed a peak— 
between-means counting technique in which the exceedances are 
accumulated from the largest value to the smallest. The two 
methods are illustrated in figure 7. Appendix A, figure A-4, shows 
a comparison of the two counting methods on an acceleration time 
history obtained from a NASA B 57 test flight. These results 
indicate that the two techniques give the same values at the end 
points of the load factor range, and are within a factor of 2 in 
the midrange. This same behavior is also shown in the results on 
page 32 of reference 10. 
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data 


Gust and Maneuver Acceleration Separation 
Inspection of many power spectra of the eg normal acceration 
indicated that usually the low-frequency maneuver 
accelerations were sufficiently far removed from the gust responses 
so that suitable low-pass and band-pass numerical filters could be 
used to separate them out. Examples of such power spectra are 
given in figure 8. In these spectra, it is believed that the peaks 
below about 0.1 Hz are due to pilot induced maneuvers; those 
between 0.2 and 0.5 Hz are the basic airframe response to 
atmospheric turbulence in the short period mode with autopilot off. 
The significant responses between 0.7 and 1.0 Hz were identified 
as the aircraft turbulence response with autopilot on in the later 
phase of the program (Vol. II) where autopilot status was 
monitored. The 1.5 Hz response peak is believed to be due to the 
wing first bending mode. Accordingly, the filters illustrated in 
figure 9 and described in Appendix D were developed based on the 
methods of reference 11 and utilized here to separate the pilot- 
induced accelerations from the aircraft gust response. The break 
frequency selected for this aircraft was 0.09 Hz; the top of the 
band-pass was set at 1.2 Hz to remove structural resonances, that 
is, wing first bending mode at 1.5 Hz (fig- 8). Results of the 
application of these filters to a typical time history are given 
in figure 10. 

In the earlier analog VGH program, the gust peaks were 
identified by eye by their sharp rise times compared to the 
maneuver- induced g loads, and the increment of the gust 
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acceleration relative to the maneuver acceleration was measured 
directly from the time history trace (ref. 12) . Due to the large 
amounts of manual labor involved, the minimum acceleration 
increment utilized was usually 0.2 g, and occasionally 0.4 g 
whereas in the present DVGH program, g increments of 0.05 g have 
been utilized at small absolute g levels. 


Gust Velocities 

Derived gust velocities U de were computed using the method 
of reference 13, and the band-pass component of the normal 
acceleration at the eg, a nG . Thus, 


U 6e= 


9W 

K ,p. C L SV. 


nG 


where K g = 


0.88 

5.3 + fi 


Mg 


and M g = 


2W 


C L p cgS 


In the current program, the lift-curve slope C L is the 

a 

untrimmed flexible lift-curve slope for the entire airplane, and 
is a function of Mach number, altitude, and flap deflection and 
is given in Table II. Time histories of were computed using 

the above equation; U de exceedances were then determined using the 
level crossing technique on these time histories. 


Statistical Formats 

Much of the data are given as a percentage of the total flight 
time that is spent in some particular condition, or combination of 
conditions. Thus, altitude is broken down into 9 bands (or bins 
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or intervals) 5000-feet thick, and the time spent therein is 

reported as a percentage of the total flight time (201 hours) . 

(Some sorts were done using data base of 201.4 hours, some using 

a base of 201.13 hours.) The bands for the major parameters are: 

Pressure altitude 5000 feet 

Airspeed 1° knots 

Mass 30,000 pounds 

Duration 0.5 hours 

Flap deflection detent 

In other instances, data are given as a percentage of total 
flights for a particular condition such as trip length, maximum 
altitude, and maximum normal acceleration. 

RESULTS AND DISCUSSION 
Organization 

The detailed numerical results are presented in two sections: 

1. Flight Profile Statistics 

2. Acceleration Derived Statistics 

Flight Profile Statistics are given on a Percent of Total 
Flight Time (201 hours) and Percent of Total Flights (83) basis. 
Some Flight Profile Statistics are given for the entire flight 
(flaps up or down) , for flaps down only, and for spoilers deployed. 
The Acceleration Derived Statistics are given on a Counts/Hour and 
Percent of Flights basis. The groupings for both results are given 
in figure 11 together with their corresponding figures 12 through 
24. The differences between the DVGH and VGH systems are 
summarized on figure 25, and figure 26 attempts to show the effects 
of several of these differences on the exceedance results. 
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In addition, some effects of autopilot operation on normal 
accelerations are discussed in the section "Autopilot Effects" and 
illustrated in figure 27. 

Flight Profile Statistics 

Many Flight Profile Statistics hitherto not generally 
available have been compiled using a subset of the data types 
available on the original DFDR. These statistics are discussed 
here for the three conditions shown in figure 11: (1) for the 
"Entire Flight" from takeoff to landing with flaps up or down, 
(2) for "Flaps Down Only," and (3) for "Spoiler Deflections." 

Entire Flight .- After several iterations with airline 
aircraft designers, it was decided to provide the Flight Profile 
Statistics for entire flights in the following ways: 


12(a) 


Gross weight histograms for takeoff and landing 


12(b) 


Fuel weight at takeoff and landing versus trip 
duration: matrices 


13(a) , (b) 
(c) 


Gross weight and altitude times for flight modes of 
climb, level, and descent: matrices and plots 


14(a) , (b) 
(c) 


Airspeed and altitude times for flight modes of 
climb, level, and descent: matrices 


15(a) , (b) 
(c) 


Maximum altitude per flight versus flight duration: 
matrix and plots 


The plots do not include all data in the parent matrices; 
rather they represent summary trends of interest. The matrix 
formats themselves show all available data and are constructed 
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to permit rapid identification of areas of maximum activity by a 
visual scan without plotting. 

Figure 12(a) shows that no flights took off weighing more than 

430.000 pounds; over 40 percent weighed between 370,000 to 400,000 
pounds. The smallest takeoff mass, 280,000 to 310,000 pounds was 
used for about 7 percent of the flights. In landing, none of the 
aircraft grossed less than 280,000 pounds, nor more than 370,000 
pounds. The matrices of fuel weight at takeoff and landing versus 
trip duration, figure 12(b), show the most popular trip length 
(31.3 percent of the flights) was 2 to 2-1/2 hours carrying 70,000 
to 100,000 pounds of fuel at takeoff; 22.9 percent of flights 
landed with between 10,000 to 40,000 pounds of fuel. For all trip 
lengths, 57.7 percent of the flights landed with between 40,000 to 

70.000 pounds of fuel. 

The joint distributions of total flight time spent in weight 
and altitude bands for climb, level, and descent flight modes are 
given in matrix form in figure 13(a) and plotted in summary form 
n figure 13(b). The matrix shows that in climb, 6.5 percent of 
the time is spent at 370,000 to 400,000 pounds with each 5,000 feet 
riltitude band up to 34,500 feet showing roughly the same 1 percent 
usage. In level flight, the most frequent condition is between 

340.000 to 370,000 pounds gross weight at 33 percent of the time. 
At these conditions, 17 percent of the total time is spent between 
29,500 and 34,500 feet, with the next higher and lower altitude 
band carrying 9 percent and 5 percent, respectively. In descent, 
7 percent of the time is at 310,000 to 340,000 pounds, with each 
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of the three altitude bands below 14,500 feet showing 1 to 1.7 
percent of the total flight time. The summary column on the right, 
plotted in figure 13(a), shows that about 16 percent of the time 
is spent in climb at all weights and altitudes, 69 percent in level 
flight, and 15 percent spent in descent. The summary rows for time 
spent in each altitude band for climb, level, and descent flight 
modes are plotted in figure 13(c). 

Joint distributions of the total flight time spent in 10- 
knot CAS intervals and 5000-feet altitude bands are given in matrix 
form in figures 14(a), 14(b), and 14(c) for climb, level, and 
descent flight modes, respectively. Figure 14(a) shows that in 
climb, CAS as high as 370 to 380 knots were used in the altitude 
range 9500 through 24,500 feet for 0.0093 percent of the total 
flight time (about 1 minute). And for 3.2 percent of the total 
flight times, climb airspeeds of 340 to 360 knots were used in 
traversing this altitude band as indicated by the dotted boxes. 
The lower CAS used above and below this altitude band can be 
readily identified as shown in the remaining dotted boxes in the 
figure. In level flight, the most prevalent conditions are 
CAS = 300 to 320 knots at 29,500 to 34,500 feet about 25 percent 
of the time, and 290 to 300 knots at 34,500 to 39,500 feet about 
12.5 percent of the time. Descent shows a broader distribution; 
the single most popular condition is 140 to 150 knots in the lowest 
altitude layer in terminal cruise, approach, and landing condi- 
tions. These matrix-type plots readily lend themselves to hand- 
sketching-in of contours of constant percent time; an example of 
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this is given in figure 14(c), for descent, where a contour 
enclosing all points greater than 0.25 percent of time has been 
indicated by dashed lines to illustrate the broad distribution for 

this parameter. 

Figure 15(a) gives the matrix of percent of flights to maximum 
altitude versus flight duration. Plots in figures 15(b) and 15(c) 
show that more than 75 percent of all flights went to a maximum 
altitude greater than 29,500 feet. Most of these had flight 
durations greater than 2 hours. About 23 percent of the flights 
went to a maximum altitude between 34,500 and 39,500 feet for trips 
2 to 2-1/2 hours long. About 5 percent of the flights had maximum 
altitudes less than 9,500 feet and durations of one hour or less. 

Flaps Down Only .- Flap detent position data for the trailmg- 
edge flap surface are given in figure 16. Note that any flap 
deflection within the detent limits shown in figure 16 were 
categorized as "in the detent." On takeoff, the initial setting 
was 10 -degrees detent (except in two cases which took off with 22 
degrees) lasting about 0.13 percent of total time, followed by 4 
degrees for about 0.70 percent of total time. By definition, the 
takeoff phase begins at lift-off and ends the first time the flap 
setting goes to zero. Subsequent operations with flaps down are 
in the landing phase. Flap deflection above 10,000 feet did not 
usually occur, except for one instance where 3-degree flap was 
selected while in a holding pattern at 15,000 feet at 210 knots 
CAS. Fourteen minutes later, the aircraft descended through 10,000 
feet and the data were picked up in the regular computer analysis. 
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For landing, the most used detents were 4, 10, 22, and 42 degrees; 
the 18, 27, and 33 degree detents were generally transited rapidly. 

For each of the most used detents, the gross weight, altitude 
above airport, and airspeed distributions are given separately for 
takeoff and landing, figures 17(a) through 17(f). For reference, 
the flap placard limit speed is also shown on the airspeed 
distributions. It can be seen that the airspeed distributions 
become more sharply peaked as the flap deflection increases. 

Spoiler Deflections .- Data showing the operation of spoilers 
2R and 5L (see fig. 1 for locations) are given in matrix form as 
a percent of total flight time in a given deflection band and 
within a given airspeed band in figures 18(a) and 18(b) and in 
plotted form in 18(c). Altitudes above which spoiler deflections 
are greater than 10 degrees are plotted in figure 19. These data 
indicate that most spoiler operations occur at about 240 to 260 
knots at altitudes between 4,500 to 14,500 feet. They also show 
that spoiler 2R is used about 10 percent more often than spoiler 
5L. These usages are as speed brakes only, with flaps up. Note 
that spoiler 2 moves linearly with speed brake handle position. 
Spoiler 5 moves nonlinearly with handle position, such that the 
spoiler angle for number 5 is much less than for number 2, until 
maximum spoiler deflection angle is approached. In addition, 
deflection of individual spoilers is limited by available control 
hinge movements; the full spoiler deflection of 60 degrees normally 
can be reached only at airspeeds below about 200 knots. In the 
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present data set, the maximum deflection of 55 to 60 degrees was 
reached only briefly (3 to 4 seconds at 200 to 220 knots). 


Acceleration Derived Statistics 
The acceleration level crossing counts per hour results, 
obtained within the given pressure altitude bands, are given as 

follows: 

Quantity Data Matrix 


Plots for altitude 


nM 

i rvG 


fig. 20(a) 
fig. 20(b) 
fig. 20(c) 
fig. 21(a) 


j-d,e,f,g,h,i,j,k,l 

b,c,d,e, D 


The following observations are noted for the normal acceleration 
results in figure 20: 


1. 


The a n = 
altitude. 


0 level crossing counts/hour increase with 
generally . 


2 . 


The a = 0 level crossing counts/hour for a nG are 
slightly higher than those for a n because of the 
biasing induced in the a n by the positive load factor 
in maneuvers as shown in the sketch below. 


Raw Data 


Maneuver Only 


Ar 


a n 0 


-.4 


jrA. 




_L 


T i nu^ 


.4 

a nH 0 

-.4 


J£=L L 


_L 


T 1 TDP 


.4 

a nG 0 

-.4 


Gust Only 





J L 


T i mo 


3. Approximate positive-negative symmetry of the data is 
observed 

The a y data, figures 21(a) through 21(j), show slightly 
higher rates at a y = 0 at the higher altitudes. 

The U de exceedances derived from the a nG are given in 

figures 22(a) through 22 (j). 
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The maximum positive and negative normal accelerations 
experienced on each flight are shown in figures 23(a) through 
23(f). The results indicate that for a n , about 30 percent of the 
flights experienced between 0.25 and 0.30 plus g's, and 37 percent 
experienced between -0.20 and -0.25 g’s. The most frequently 
experienced a y was in the .04 to .06 interval both positive and 
negative. The 0.15 to 0.20 maximum gust g's per flight, 
figure 23(d), appear reasonably symmetric at about 35 percent of 
the flights. Interestingly enough, the most prevalent negative 
maneuver g's, figure 23(e), occurred significantly more often 
(57 percent of the number of flights) than the most prevalent 
positive (39 percent of the number of flights) . These maximum g 
level data have been compiled in figure 23(f) to show the percent 
of flights to exceed a given g level. Thus, +0.3 a n was 
exceeded on about 45 percent of the flights, +0.3 a^ on about 5 
percent of the flights, and +0.3 a nG on 17 percent of the flights. 

For each flap detent position, the a n level crossing counts 
per hour are given in matrix form in figures 24(a) and 24(b) for 
take off and landing. Corresponding plots are given in figures 
24(c) and 24(d). The dotted line in those plots is taken from 
figure 20(d) for a n . The only significant difference in these 
results is that for takeoff with 4 and 10 degrees of flaps, 
negative normal accelerations are experienced somewhat more 
frequently than for the data from figure 20(d), which are for any 
flap setting, and at +0.4g, where the data indicate higher rates 
than with flaps down 4 and 10 degrees. 
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Comparison of DVGH and VGH Data 

Previous published information by NACA/NASA on flight load 
experiences of transport airplanes during routine operations (i.e. 
refs. 1-5) were derived from analog type velocity-acceleration- 
height (VGH) recorders using manual data reduction techniques. The 
data in this report were obtained from digital recorders utilizing 
an automatic data reduction process. In both cases, the primary 
measurement is acceleration measured very near the airplane's 
center of gravity. Some of the differences between the two 
recording and processing techniques are listed in figure 25. 
Before comparing results, a short synopsis of the two systems is 
given. 

VGH Program .- The frequency response of the VGH recorders were 
from near DC to 5 Hz. Reading accuracy and overall system errors 
combined to provide an accepted band width of 0.01 to 5 Hz for gust 
and maneuver load determination (ref. 14) . Data were recorded as 
a continuous acceleration time history on 70mm photographic paper 
roll film. Peak accelerations were manually read in +0.01 g 
increments from the 1.0 g level flight reference line. Only the 
maximum peak occurring between successive crossings of the 1.0 g 
Level flight reference line were read and counted (see fig. 7) 
Small oscillations of magnitude +0.05 g or less (and occasionally 
up to +0.3 g) about the zero level were not counted. The peaks 
were then cumulatively totaled from the highest level to the lowest 
to produce a frequency of exceedance distribution. Separation of 
gust and maneuver accelerations were dependent on the film reader's 
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experience but, generally, gusts were identified by their rapid 
rise and decay time, i.e., they had a relatively higher frequency 
than maneuver loads. The method of manual reading and 
representation of the analog data became known as the cumulative 
peak-between-means counting technique. 

Digital Flight Data Program .- The process for reducing data 
from the digital recorders was more automated and, therefore, 
differs somewhat from that in the VGH Program. First of all, the 
sample rate was limited to 4 per second since that was already in 
use on the DFDR system. The maximum usable data frequency, 
therefore , would be near 2 Hz. Appendix A addresses the effect of 
sampling rate on determination of turbulence induced accelerations. 
It shows that during a 4.5-minute turbulence encounter of a B 57 
airplane equipped with a high response gust measurement system 
(20 Hz accelerometer cutoff) there was a loss of 1/3 to 1/2 of the 
acceleration peaks as the data reduction sampling rate dropped from 
40 per second to 20 per second to 4 per second. Also, the DFDR 
acceleration data were filtered to reject frequencies above 1.2 Hz 
to eliminate wing— body elastic response. Automatic separation of 
gusts and maneuver accelerations were accomplished using a low pass 
(maneuver) filter, DC to 0.09 Hz, and a high pass (gust) filter, 
0.09 to 1.2 Hz. 

Comparison of Results .- Figure 26 compares the derived gust 
velocity experience from 200 hours of L 1011 data with the VGH data 
from reference 4. Also shown are data from Volumes II, IH» and 
IV of this report that included additional L 1011 data, B 727 and 
B 747 data. The two different L 1011 aircraft agree very well 
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despite an almost factor of 10 difference in the number of flight 
hours. The L 1011 data are noticeably higher at the lower gust 
values than the Boeing aircraft and the VGH gust curves derived 
from reference 4. All of the curves seem to converge at about 
20 feet per second gust velocity except for the extrapolated B 747 
curve from volume IV. The consistently lower values for that 
particular aircraft may be explained by the fact that it was used 
almost entirely on long overwater flights during the period of data 

collection . 

The higher slope of the digital data may be related to 
differences between the manual and automatic data reduction 
processes. In the manual reading of VGH data, the acceleration 
peaks are increasingly more difficult to distinguish as they become 
smaller in the photographic time history traces. Thus, there may 
be a tendency to undercount VGH accelerations, whereas the digital 
process would provide an accurate count regardless of the peak 
magnitude, provided the sample rate were adequate. As previously 
discussed, the 4 sample per second rate used to obtain the present 
digital results may have missed counting a significant number of 
jC -ual acceleration peaks. If this is the case, then the digital 
xceedance counts would be higher than shown. 

Because of the many differences between the two data 
collection programs and because the current digital data suggests 
a change in slope of the derived gust velocity experience from 
previous VGH results, the author does not feel that it is practical 
to define a correlation factor to permit combining all of the data 
sets. It is recommended, however, that any future digital flight 
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loads program for transport airplanes consider increasing the 
frequency response of the gust measurement and analysis system to 
a level that, as a minimum, matches the VGH system. The new system 
should provide information at not only the aircraft short-period 
response frequency, but also at the higher frequencies associated 
with significant aeroelastic modes of the aircraft structure. 

Autopilot Effects 

Although autopilot status was not monitored in this test, some 
confirmation of some of the effects predicted by reference 15 was 
gleaned from the present results. Thus, the theoretical results 
of reference 15 indicate that the effects of the autopilot should 
be: 

1. "...the introduction of multiple response modes at 
frequencies both below and above the controls-f ixed 
short-period frequency. . . , " and, 

2. to reduce the gust response magnitude by 10 to 
25 percent. 

Some evidence of the existence of a resonance at a frequency 
below the controls — fixed short-period frequency may be seen in the 
sporadic appearance of a low amplitude limit cycle oscillation on 
a . A typical occurrence is shown in figure 27. The occurrence, 
magnitude, and frequency of this oscillation are more fully 
documented in the test reported in Volume II where autopilot status 
was monitored. The effect of the autopilot on the response above 
the controls-f ixed short-period frequency is shown in figure 8 
where the data intervals evidently overlapped autopilot "off" and 
"on" periods, resulting in some response at -0.2 Hz and some at 
-0.8 Hz. This too is more fully illustrated in Volume II. 
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However, the effects of the autopilot operation on the magnitude 
of the gust response and its effect on a n and the derivation 
could not be determined here, due to the lack of autopilot status 
and other information. Autopilot effects were not accounted for 
in the results of the VGH program, references 1-4. 


CONCLUSIONS 

Initial results of an effort to utilize a limited data set 
from existing aircraft Digital Flight Data Recorders to describe 
the aircraft operating conditions indicate: 


1. A significant problem in Digital Flight Data Recorder 
data processing is data editing, particularly identifica- 
tion and replacement of wild points. A two-parameter 
algorithm has been developed and successfully applied to 
this problem, replacing the manual methods initially 
developed . 

2. A large variety of Flight Profile Statistical Data useful 
to airline aircraft mission analysts and designers can 
be compiled from the few parameters selected. 

3. An objective technique has been developed and applied to 
separate the maneuver and gust components of the normal 
acceleration data. 

4. Acceleration gust exceedances derived from the DFDR 
system at 4 samples per second may be significantly less 
(approx. 50%) than if actual peak values were counted. 


5. Improved data systems with sample rates and accelerometer 
natural frequencies higher than those used in the DFDR 
system will be needed to adequately describe the aircraft 
response to atmospheric turbulence. 
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APPENDIX A 


ACCURACY CONSIDERATIONS 

A complete study of the accuracy of each data type was not 
undertaken. However, due to its importance, some checks were 
performed on normal acceleration as described below. These results 
are also generally applicable to the lateral acceleration results 
given herein. 

Accuracy 

A basic DFDR system (ref. 6) was flown on a NASA DeHaviland 
Twin Otter aircraft to assess overall system performance by 
comparing its normal acceleration output with that of a NASA 
measuring and recording system. The overall accuracy of the NASA 
system was +0.12 g's, or about twice as accurate as the DFDR system 
at +.2 g's. The results of the test, shown in figure A-l, show 
that the DFDR system results had a standard deviation of about 
.055 g relative to the NASA system, with a correlation of 0.9790. 
In addition, figure A-2 indicates that the power spectrum from the 
DFDR system agrees well with that of the NASA system, with the 
largest discrepancy at around 0.02 Hz. Although data were obtained 
in a short time on only one flight, they indicate that the basic 
DFDR normal acceleration data are being accurately measured and 
recorded. 

Sample Rate 

The DFDR system described in reference 6 usually provides 
normal acceleration data at 4 samples per second. This provides 
adequate frequency response to about 2 Hz (based on the Nyquist 
rule) , and was judged adequate to define the principal rigid-body 
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Normal Acceleration, g s 
NASA System 


.Figure A- 


a = .98865665 a + -00611^16>752 
n DFDR n NASA 


Standard deviation = +.054772 g's 

> 

Correlation Coefficient = .9790 


Comparison of DFDR system normal acceleration measurements 
with NASA system results. 
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FIGURE A-2 . - Comparison of normal acceleration power spectra 
measured by DFDR system and NASA measuring and 
recording system on a NASA DeHaviland Twin Otter 
aircraft. 



gust response of large transports and for computing derived gust 
velocity, U^. (Some data were obtained at 8 samples per second. 
The spectra for 4 and 8 samples per second, shown in figure A-3, 
are in good agreement up through about 2 Hz.) The question arises, 
however, about the effect of digital sampling on peak count since 
for previous analog VGH data the actual peaks were read — in other 
words — what is the reduction in peak count due to digital sampling? 
A recent unpublished analysis of a 200 samples per second digital 
eg normal acceleration record from a 1982 flight of a NASA B 57B 
aircraft in moderate turbulence is summarized in figure A-4 . The 
results show that exceedances (both level crossing and peak count) 
increase significantly with sample rate up to 20 samples per second 
and only a slight increase when sample rate is increased from 20 
to 40 samples per second. These results would indicate that 
exceedances reported herein may be 1/2 to 1/3 of those that would 
be determined from the previous method. This should be considered 
in the application of the acceleration data presented in Volumes 
l through V of this report. 

If higher frequency response data are important in the future, 
consideration should be given to accelerometer frequency response 
haracteristics. The so-called cutoff frequency (3 dB attenuation) 
typical to DFDR usage is 4 Hz. The accelerometer utilized for the 
research program with the B 57B aircraft had a cutoff frequency of 
10 Hz. 
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Figure A-4.- Comparison of acceleration exceedances obtained from unpublished NASA B57 4.5 minute run number 
3113 using Level Crossing counting technique at 4, 10, 20, and 40 samples per second with exceedances 
obtained using the cummulative Peak-Between-Means counting technique. 



APPENDIX B 


WILD POINT EDITING ALGORITHM: KEDIT 

This algorithm was developed in 1979 by D. A. Keskar of System 
Development Corporation Integrated Services, Inc. as part of 
contract NAS1-15400. The description given herein is based on that 
Company's document PDD-79-01, titled Program Description Document 
for Automatic Edit Program for Digital VGH Data Analysis. 


OBJECTIVE 

A software package needed to be developed for automatic data 
editing. It had to conform to the following restraints. 

1. It should not replace any "good" data, with no artificial 
limit imposed on the magnitude of data excursions. 

2. it should replace a "bad” point with the most probable 
data value at that point. 

3. It should reconstruct data gaps caused by loss of frame 
synch. 

4. The software package must be optimized as it will be used 
in editing hundreds of tapes, each 25 hours long, with 
sample rates as high as 4 per second. 
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PROBLEM SOLVING 

The idea used in solving the problem is from basic statistics. 
Since the data are collected from the aircraft flying in a 
realistic physical environment, it must meet the tests based on 
moving average and standard deviation (fig. B-l) . 

The algorithm used is as follows: 


For 

calculate 

and 


i = p+1, P+2, 

P 


n-p 


x < = 

2p 


y *h, 

Lt D=-P 

(j^O) 


.2 


v^ +p - 

-L. t >, x i 

2p 3=~ 


i+j 

P 


] - ( X, ) 2 


mean 


variance 


then 

e i = x i " *1 

If | e | > K • Sj perform "local mean" test 

If | e . | < k • S i Xj is left unchanged. 

The key to the success of this algorithm lies in proper selection 
of values of p and K. 

Certain modifications must be done in the algorithm before it 
can be implemented. It should be noted that the concept of moving 
average fails if the point to be edited occurs just before a data 
drop out. To eliminate this problem, the points are temporarily 
replaced by the average of the previous 10 points. This results 
in a continuous record of data to be edited, (fig. B-2) . 
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An analytical investigation was undertaken to find the most 
suitable values of constants p and K . As shown in figures 
B-3(a) and B-3(b), the decrease in the value of p below some 
threshold results in improper editing as the statistical informa- 
tion available is insufficient; while an increase in the value of 
p results in additional computational burden, as is obvious from 
the algorithm. For a reasonably good record (fig. B-4) , values of 
p > 4 is sufficient as shown in figure B-3(a). On the other hand, 
for a record with a large number of noise spikes and out-of-synch 
points, figure B-5, p should be at least 8 as shown in figure B- 
3(b). Since the algorithm is used for editing many tapes, each 25 
hours long, computational time requirement was a major 
consideration. For the data records considered, any value of p 
greater than 8 will edit the data satisfactorily, but will add 
significantly to computational cost by increasing the execution 
time. Based on this study, a value of p = 8 was found most 
suitable. 

Similarly, a small value of K will replace some good points 
in the data while a large value of K will pass noise spikes in 
the data. The results of varying K in editing a typical data 
record are shown in figure B-5(a), B-5(b), and B-5(c). Thus, 
a few noise spikes were passed as good data with K = 4, while K 
= 3 replaced all the "bad" points by the most probable value at 
that point. This analysis resulted in choosing K = 3, which will 
insure replacing all "bad" data points while leaving most "good" 
data points unchanged. 
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The final algorithm computes mean and variance based on 8 
points before and after the point in question. Since the points 
"after" are not edited as yet, the presence of noise spikes will 
add significant bias in mean and variance. To take care of this 
situation, the program checks the "after" points for unusually high 
values; for a point to be included in the statistics computation, 
it must be less than 25 times the value of previous point. 

The objective of this algorithm is to replace all the bad 
points and retain almost all "good" data points. In a situation 
where mean and/or variance is close to zero, some "good" points 
may fail the edit test and eventually will get replaced. To 
alleviate this problem, one more test is done after the point in 
question fails the edit test. 


Compute : 


X i-1 I + 


M -2 


(local mean) 


If | Xj | >2 ' x L ? replace Xj with x 1 
If | Xj | <2 ‘ x L ; Xj is left unchanged 


'his modification significantly improves the chance of retaining 
almost all "good" data points. 

The final flow chart for the algorithm is given in figure B-6. 
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KEDIT SUMMARY 


A software package has been developed for automatic editing 
of data obtained from airline aircraft Digital Data Flight 
Recorders. The algorithm has been tested on several data records 
to insure that all "bad" data points are replaced by the most 
probable data value at that point and almost all "good" data points 
are left unchanged. The algorithm also reconstructs the data gap 
of 4 seconds due to loss of frame synch. The algorithm is 
computationally efficient and takes about 6 seconds of CPU time to 
edit a record 10 minutes long sampled at the rate of 4 per second 
on a CDC-6000 series computer. The implementation has resulted in 
savings of 8 to 10 man hours per tape and has also improved the 
quality of editing since "human errors" are eliminated. Finally, 
the algorithm compiles statistics of types and numbers of points 
edited in each channel. This could be a vital piece of information 
in interpreting the final results. For example, if the number of 
points edited in one particular tape is too high, the confidence 
in the subsequent results obtained from the data of this tape 
should be low. 
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Figure B-1. - Illustration of algorithm. 


999999 

I 

I 

I P - Replaced Point 

I 

I 

I 

I 

I 

I 

/ 



Figure B-2. - Temporary replacement of loss of synch data- 
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(a) A "good" record (figure B-4a) ( b ) A n bad .. record (figure B-5a) 


FIGURE B-3 . - Effect of varying p on goodness of fit. 
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(b) Edited record with p=3 - no change 


FIGURE B-4 . — Illustration that algorithm does not alter good 
record. 39 
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(b) Edited record with K=4 
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(c) Edited record with K=3 


B-5 . - Illustration of effect of K on editing "bad" record 
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FIGURE B-6 


Flow chart of completed algorithm 









APPENDIX C 


WEIGHT CALCULATION 

Since the fuel weight on the L 1011 aircraft can amount to 30 

percent of the gross weight (130,000 lbs out of 430,000 lbs), it 

was deemed necessary to determine a reasonably accurate aircraft 

weight at each instant of flight to permit accurate computation of 

the derived gust velocities. To permit this, the airline operator 

supplied for each flight the following support data: 

Gross weight leaving the terminal 
Gross weight at beginning of takeoff roll 
Fuel weight at beginning of takeoff roll 
Fuel burned terminal to terminal 
Airborne time 
Taxi-in time 

It was further determined that the average fuel burn during 
the takeoff roll was 328 pounds. 

These data were used to compute the gross weight at lift-off 
and at touch-down as shown in figure C-l. The gross weight at any 
time in between was computed assuming a linear variation with time. 
This linear assumption was found to be, for one flight, within 2 
percent of the weight computed using fuel-burn equations supplied 
by the manufacturer. This comparison is shown in figure C-2. 
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Gross weight, k lbs. 


1 


Leave terminal 



Time from leaving terminal# hours 


O Airline Operator Supplied for Each Flight 
o Gross weight @ 1 

o Gross weight 0 2 

o Fuel weight 0 2 

o Fuel burned terminal to terminal = AGW t . 5 
o Airborne time = t 3 _ 4 

o Taxi-in time = t 4 _ 5 in minutes; x 100 Ibs/min - agw 4 _ 5 
o Takeoff fuel burn = 328 lbs 

0 NASA Calculated 

o Gross weight 0 3 = GW 0 2 - 328 lbs 

o Gross weight § 5 = GW ?1 “ AGW V5 

o Gross weight 0 4 = GW 0 5 + AGW 4 _ 5 

o Burn rate 3 . 4 = GW 0 3 - GW @ ..A 

t 

^ 3-4 

o Gross weight (t) = GW 0 3 - BR 3 _ 4 x ^ 


Figure C-l.- Scheme for gross weight calculation. 
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Gross weight klbs. 



Figure C-2.- Comparison of NASA staightline method of calculating 
gross weight variation with time versus Lockheed fuel flow 

algorithms methods. 




APPENDIX D 


DESIGN OF NUMERICAL FILTER TO 
SEPARATE GUST AND MANEUVER ACCELERATIONS 

The 5112 filter of reference 11 has been applied to separate 

X 

and high frequency components of the normal acceleration response 
as indicated on page 13. The application of that technique to this 
problem was performed by Dwight W. Smith, formerly of System 
Development Corporation Integrated Services, Inc., as part of 
contract NAS1-15400, and is described in that Company's document 
PDD 77-2, August 24, 1977. The design parameters and resulting 

numerical weighting functions are given in figures D-l, D-2 , and 
D-3 . 
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(a) Desired frequency response of filters and definition of 
critical frequencies. 
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(b) Actual frequencv response of filters with 62 weighting 
functions at 4 samples per second. 


Figure D-l. - Numerical filters used in separating gust and 
maneuver accelerations. 
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a nM ( fc 


a nG 

where 


Figure D-2 


i) = a n ( fc i) xW i + 

LP 


NWTS 

E" 


LP 


[an (tj.j- +1 ) + a n (t j+j ] 


j=2 


i) = a n < fc i) x w i + 
BP 


NWTS 

E" 


[a n (t i . j +1 ) + a n (t j+j .,) ] 


3=2 


tj = ith time point 

W- = filter weights W . , , given in figure D-3 
NWTS = Number of weighting constants = 62 
LP = Low Pass 
BP — Band Pass 


- Time history weighting function format. 
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Wj — w, , w 2 r W 3 , ... 

w « 

LP 


1 

. 49463E-01 

16 

. 13630E-01 

31 

- . 79834E— 02 

47 

. 29154E-02 

2 

. 49244E-01 

17 

. 10456E-01 

32 

— . 74532E— 02 

48 

. 30151E-02 

3 

. 48588E-01 

18 

. 74454E-02 

33 

-. 67819E-02 

49 

. 30230E-02 

4 

. 47509E-01 

19 

. 4 63 57E-02 

34 

-. 60000E-02 

50 

. 29482E-02 

5 

. 46023E-01 

20 

. 20579E— 02 

35 

-. 51380E-02 

51 

. 28011E— 02 

6 

. 44 158E-01 

21 

-. 26264E-03 

36 

— . 42255E-02 

52 

. 25933E-02 

7 

. 41943E-01 

2 2 

23065E-02 

37 

- . 32909E— 02 

53 

. 23367E-02 

8 

. 39418E-01 

23 

— . 40603E— 02 

38 

- . 23606E— 02 

54 

. 20440E— 02 

9 

. 36623E— 01 

24 

-. 55171E-02 

39 

-. 14584E-02 

55 

. 17274E-02 

10 

. 33 605E-01 

25 

- . 667 58E-02 

40 

- . 60535E-03 

56 

. 13985E-02 

11 

. 304 13E-01 

26 

- . 754 10E-02 

41 

. 18061E— 03 

57 

. 10685E— 02 

12 
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27 

- . 81231E-02 

42 
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58 

. 74719E-03 

13 

♦ 23711E-01 

28 
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43 

. 14967E-02 

59 

.44321E-03 

14 

. 20305E-01 

29 

- . 8 5027E-02 

44 

. 20083E-02 

60 

. 16377E-03 

15 

. 16928E-01 

30 

-. 83427E-02 

45 

. 24154E— 02 

61 

- . 85417E-04 
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. 27171E-02 

62 

-. 30021E-03 


(a) Low Pass 


W l 

BP 


1 
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16 

-. 23620E-02 

31 

. 82254E— 02 

47 

-. 22518E-03 

2 

.26116 

17 
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32 

. 01527E— 01 

48 

. 44415E— 02 

3 
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18 
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33 
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49 

. 34913E-03 
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21 
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36 

. 26061E— 02 

52 
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7 

-. 84234E-01 

22 
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37 
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53 
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8 
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23 
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38 
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54 

52614E-03 

9 
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24 
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39 
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55 
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10 
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25 
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40 

. 43064E-05 

56 
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11 
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26 

- . 4 1775E-02 

41 
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57 

- . 22346E-02 

12 

- . 28227E-02 

27 

- . 65077E-02 

42 

. 67874E— 02 

58 

. 3 1915E— 03 

13 

- . 44010E-01 

28 

. 104 02E-01 

43 

- . 59443E-03 

59 

- . 20560E-02 

14 

-. 33080E-01 

29 

. 48957E-02 

44 

. 32631E— 02 

60 

- . 29623E-02 

15 

. 13860E-03 

30 

- . 45374E-02 

45 

. 63120E-02 

61 

- . 25083E-03 
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- . 28507E— 03 

62 

-. 10833E-02 


(b) Band Pass 


Figure D-3 Filter weights 
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TABLE I.- WEIGHT AND GEOMETRY OF THE L 1011 AIRCRAFT. 


Weights 

Maximum takeoff 
Empty weight 


430,000 lbs 
247,500 lbs 


Areas 


Wing 3,456 ft 2 

2 

Stabilizer and elevator 1,282 ft 

Mean chord 

Wing 22.3 ft 

Stabilizer and elevator 19.42 ft 

Sweepback quarter chord 

Wing 35 ° 


Stabilizer and elevator 

Fuselage Stations of Mean Aerodynamic Chord Leading Edge 


Wing 1143 

Stabilizer and elevator 1885 

Accelerometer location (see figure 1) 

Fuselage station 1243 

Waterline 182 
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TABLE II.- LIFT-CURVE SLOPES USED IN CALCULATING U de 



FROM FOR 

THE L-1011- 

1 AIRCRAFT 



LIFT- 

-CURVE SLOPE 

C L , PER 

a 

DEGREE 


Flaps 

up C L = f(M 

d 

t HP) 


M 

HP = 0 

10 

20 

40 kft 

.20 

.0923 

.0928 

.0929 

.0936 

.35 

.0923 

.0928 

.0930 

.0938 

.50 

.0913 

.0920 

.0929 

.0946 

.60 

.0918 

.0928 

.0940 

.0963 

.70 

.0940 

.0954 

.0970 

.1003 

o 

00 

• 

- 

. 1038 

.1058 

.1100 

.89 

- 

.1210 

.1240 

. 1305 

.91 

- 

- 

.1227 

.1286 

.95 

- 

- 

.1030 

.1081 


Flaps 

Down C L = f(FLP) 



or 


FLP.dea HP = 0 


0 

.0925 

4 

.0973 

10 

.0980 

18 

.0975 

22 

.0971 

27 

.0962 

33 

.0948 

45 

.0912 
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FIGURE 1.- Aircraft three-view with locations of accelerometers and spoilers. 
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Figure 2. - Location of service area and scope of data base. 




AIRLINE READOUT AND TRANSCRIPTION: APPROXIMATELY TWICE PER WEEK 
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LANDING WEIGHT 
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Figure 3. - Digital VGH program data sources and handling. 
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Figure 4. 
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Simplified data flow 
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Note that the zero level incremental 
(= l.Og absolute) is only counted 
once, that is when crossed with a 
positive slope. Also, these data 
have had all frequencies above 
1.2 Hz removed, so that zero level 
j itter is reduced without using a 
threshold. 


In the previous analog VGH program, 
threshold values of +0.2 and +0.3 g 
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Figure 7. - Level crossing and peak counting techniques. 
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Figure 8* - Sample L1011 normal acceleration power spectra. 





Figure 9.- Frequency response of numerical filters used to separate gust and maneuver 
accelerations, and to eliminate elastic responses above 1.2 Hz. 
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Figure 10.- Filter separation of normal acceleration time history. 



• Flight Profile statistics 


• Percent of Total Flight Time 


Percent of Total Flights 






o Entire Flight 
(flaps up or down) 

o Fig. 13- GW & HP for 
Flight Modes 
o Fig. 14- CAS & HP for 
Flight Modes 

o Fig. 12- Weight Stats 
§ TO & landing 
o Fig. 15- Max HP vs. 

duration 

o Flaps Down Only 

o Fig. 16- Flap Detents in 

take off and landing 
o Fig. 17- GW, AGL, CAS vs. 
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o Spoiler Deflections 

o Fig. 18- Spoiler deflections 
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Figure 11- organization of the Numerical Results and Discussion. 



Porcont of flight* 


Gross weight at takeoff rotation: 83 flights 



Gross weight at landing touchdown: 83 flights 



(a) Percent of take-offs and landings made at various gross weights 


Figure 12.- Weight statistics for take-off and landing. 
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(b) Percent of flights at various fuel weights and durations for takeoff and landing 


Figure 12. - Concluded. 
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Figure 13. - Gross weight and altitude statistics for climb, level, and descent Flight Modes 
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Figure 13.- Continued. 
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Figure 13.- Concluded. 
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Figure 14. - Airspeed and altitude statistics for climb, level, and descent Flight Modes. 
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Figure IS. - Maximum altitude and flight duration statistics. 
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Figure 17.- Continued. 
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(c) Landing with flaps = 4 deg. (1.4914 hours) 
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Figure 18. - Percent time spoiler is deflected vs. airspeed 
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(b) Level crossing counts per hour within pressure altitude bands 


Figure 20. - Continued 
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(a) By Level crossing counts per hour within pressure altitude bands 


Figure 21. - Lateral acceleration exceedances. 
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Element 

o Measurement & Recording 
System 

o Filtering 

o Maneuver & Gust 
Separation 

o Gust Counting 

o Weights 

o Lift-Curve Slope 

o Autopilot Effects 
o Total Hours 


Analog 

o Continuous with 5 Hz 
accelerometer cut-off 

o None 

o By eye 

o Peaks by eye 

o Maximum weights or 
rough estimate 

o Compressible rigid 
wing-only theory 

o Not monitored 

o 25,000 (ref. 4) 


Digital 

o 4 samples per second with 
4 Hz accelerometer cut-off 


o Filter cut-off at 1.2 Hz 
(to eliminate elastic 
response) 

o By computer, defined by 
frequency 

o Level crossings by computer 


o Calculated (believed to be 
within 2%) 


o Compressible elastic total 
airplane provided by 
manufacturer 


o Monitored but not accounted 
for 


o 5,000 


Figure 25.- Comparison of analog VGH and digital VGH systems. 
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Figure 26-. Comparison of results 
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Figure 27. - Low frequency oscillation on normal acceleration. 
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